Shape transitions in proton-neutron systems 
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Abstract 

We explore the possibility of shape transitions in proton-neutron systems 
driven by deformation differences between proton and neutron fluids. Within 
the framework of the proton-neutron interacting boson model, we show that 
such dynamic shape transitions cannot occur in well deformed nuclei but are 
a possibility in transitional nuclei. Likely candidates in the Os-Pt isotopes are 
discussed and predictions of the model are compared with the existing data. 
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Studies of shape transitions in nuclei at low excitation energies have been mostly 
based on coexistence of different shell model configurations . In collective models, 
shape transitions could also arise from the competition between two macroscopic 
shape variables. An example of such a shape transition, arising from the competi- 
tion between the quadrupole and hexadecapole shapes, was given in the sdg boson 
model earlier 0. Two- fluid models, where proton and neutron degrees of freedom 
are treated separately, could furnish another example. Measurements of neutron 
deformation using pion beams indicate that there could be sizable differences 
between the proton and neutron deformations. There is also indirect evidence from 
g-factor measurements that the quadrupole operators for protons and neutrons 
could be rather different, especially in transitional nuclei [^. In this letter, we in- 
vestigate in the framework of the proton-neutron interacting boson model (IBM-2) 

1^ whether such differences could also drive a shape transition in proton-neutron 
systems. For this purpose we use the analytic expansion method and the ex- 
act numerical diagonalization results from the code NPBOS 0. The expansion 
is based on angular momentum projected mean field theory and therefore provides 
an intuitive picture for shape transitions. Further, due to their analytic formula- 
tion, a global study of the parameter space can be easily carried out, which shows 
that transitional nuclei provide the most likely candidates for such a dynamic shape 
transition. As the expansion results are not very reliable for transitional nuclei, 
we use exact diagonalization in exhibiting the shape transition and applications in 
this region. 

We employ the simplest IBM-2 Hamiltonian suggested by microscopies [^, 

H = e^Ud^ + e^Udu - kQ^ ■ + ^M, (1) 

where fidp, p = vr, i/ are the (i-boson number operators for proton and neutron bosons, 
M is the Majorana operator in Casimir form and Qp are the quadrupole operators 
given by 

Qp = [dlsp + sldp] + Xp[dldp\^^\ p = 7v,iy. (2) 
The E2 matrix elements (m.e.) are calculated using the operator 

T{E2) = e^Q^ + e,Q„ (3) 
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where e^, are effective charges, and the same quadrupole operators are used 
as in the Hamiltonian. For future reference, we introduce F-spin scalar and vector 
parameters as 

Xs = {Xtt + Xz.)/2, Xv = X7T- Xv (4) 

Although other terms are sometimes included in detailed IBM-2 studies, the Hamil- 
tonian (|I|) is found to give an adequate description of level energies and electromag- 
netic (E2 and Ml) transitions |^. 

The expansion solutions follow from using the boson condensate 

\N^,K) = (iVJiVJ)-V2(5t)A^.(5t)A^.|o), hi = {l + Pl)-'/'{sl + (3,dl,), (5) 

as a trial state in a variation after projection (VAP) calculation. Here (3.,^ and jSi, are 
the mean field deformations for the proton and neutron fluids, which are determined 
from the expectation value of the Hamiltonian (|ID 

El = {N^,N,\HP,',\N^,N,)/{N^,N,\P^,\N^,N,) (6) 

by the VAP procedure. In Eq. (|^), Pqq denotes the angular momentum projection 
operator. Note that the IBM deformation variables involve only the valence nucle- 
ons, and therefore they are much larger (/5 ~ 1 — 1.4) than the typical geometrical 
model values {(3 ~ 0.2 - 0.3) |T0|. We refer to Ref. § for details of the 1/A^ expan- 
sion method, and quote here the energy expression obtained from Eq. (P) for the 
Hamiltonian (^. To the leading order in 1/A^, which is sufficient for discussion of 
systematic features, the energy surface and moment of inertia terms are given by 



El = Eo + ClL{L + 1), (7) 



(l+^2)(i"^'^^2) (8^ - 12 + (4a - 12) (x./?. + X.^) + (2a - m^X.M. 
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where Xp = ~ y 2/7xp. The quantity "a" represents the average "angular momentum 
squared" carried by a single boson and is given by 



Eq. (8) in a typical deformed (top) and transitional (bottom) nucleus. The positive 
(3 values correspond to prolate and the negative ones to oblate deformations. The 
prolate (absolute) minima are in the first quadrant and the oblate ones are in the 
third. There are 20 contour lines at 1 MeV steps, thus the maxima in the second 
quadrant are about 20 MeV high. In the deformed case, the prolate minimum in 
the first quadrant (-Emin = —4.83 MeV) is well separated from the oblate one in 
the third quadrant with an energy difference of AE = 2.23 MeV. Including the 
moment of inertia term (9) in the energy surface, this energy difference gets larger 
with increasing spin. Thus there is no chance of a cross over from the prolate to 
the oblate minimum. Although we have used a particular parametrization here for 
illustration purposes, in fact, this is a general feature of all deformed nuclei as can be 
verified by systematic studies of Eqs. (7-9). Since such a study is already available 
for IBM-1 p|, and the IBM-2 results are very similar to those of IBM-1 for deformed 
nuclei, we will not elaborate on them further. 

In the transitional case (bottom), the absolute minimum is at -Emin = —3.07 
MeV, and the energy difference between the two minima is much smaller {AE = 0.85 
MeV). Furthermore, the moment of inertia in the oblate minimum could be larger 
than that in the prolate one depending on the choice of the vector parameters in 
Eq. (4). Hence the energy difference between the two minima could decrease with 
increasing spin, with possibility of a cross over at some critical spin. Unfortunately, 
due to the soft energy surface, the 1/N expansion formulas are not very reliable 
for transitional nuclei, and a quantitative study of shape transitions is not possible 
using the analytic formulas (7-9) in this region. 

In order to exhibit the proposed proton-neutron shape transition and to study 
general conditions for its existence, we therefore rely on numerical diagonalization. 
The signature for shape change is taken from the yrast spectroscopic quadrupole 
moments, which are negative for a prolate shape and positive for an oblate one. Since 




(10) 



In Fig. 1, we show the contour plots of the energy surface Eq obtained from 
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the shape change is driven by differences in the proton and neutron deformations, 
the critical quantities to study are the F-spin vector parameters, A^v = -^tt — N^,, 
and ev, which can induce asymmetries between the proton and neutron variables. 
In principle, the Majorana interaction could also have an influence on the proton- 
neutron asymmetry. However, variations in the Majorana strength (of order 20-30%) 
did not lead to any discernible effect on shape change. In Fig. 2, we show the effect 
of these three vector parameters on yrast quadrupole moments in a typical prolate 
transitional nucleus with F-spin scalar parameters, N — 12, eg = 0.4 MeV, k = 0.15 
MeV, { = 0.17 MeV, and Xs = —0.2. The boson effective charges are taken as 
Ctt = e,^ = 0.1 eh. Fig. 2a shows the A^v dependence of the quadrupole moments 
Q{L) for Xv = 2 and ev = 0. It is seen that Q{L) changes sign between the spins 
L = 16 — 18 in the case of A^y = 4, but remains prolate in the other two cases. 
Fig. 2b shows a similar study of Xv dependence of Q{L) for A^^ = 4 and ey = 0. 
Note that the bottom results in Figs. 2a and 2b are identical due to the symmetry 
of the Hamiltonian under the simultaneous interchange of A^^ — and x-w — Xv 
Again only in the case of Xv = 2, a shape change occurs. Finally, Fig. 2c shows the 
effect of ey on the favourable case with A^y = 4 and Xv = 2. Compared to A^y and Xwi 
ey plays a marginal role. Nevertheless, it could shift the transition spin by several 
units in either direction depending on its sign; ey < makes the shape transition 
occur earlier while ey > retards it. The opposite is true when A^y and Xv are both 
negative. 

This systematic study indicates that an essential requirement for a dynamic 
prolate-oblate shape transition to occur is the coherence of the two vector param- 
eters A^y and Xv, that is they must have the same sign. The ey parameter with 
an opposite sign to the others could also facilitate the shape transition though its 
effect is less important. To understand these features better, we show in Fig. 3 
the individual proton and neutron contributions to the yrast quadrupole moments 
in the favourable (xv = 2) and the unfavourable (xv = —2) cases in Fig 2b. In 
the favourable case (top), the neutron contribution changes little from its original 
prolate value with increasing spin, while the proton contribution sharply increases, 
changing sign at spin L = 12. At low-spins, the {sd) term in the quadrupole op- 
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erator (2) dominates the m.e. and it has the same (negative) sign for protons and 
neutrons. With increasing spin, however, the expectation value of the number of 
d-bosons in the states, and hence the m.e. of the (dd) term in (2) increases while 
that of (sd) term decreases. The sign of the quadrupole moment coming from the 
{dd) part is determined by the sign of Xp which is positive for protons and negative 
for neutrons. Hence the {sd) and {dd) terms tend to cancel out for protons but add 
up for neutrons, which explains the behaviour of the individual proton and neutron 
quadrupole moments. The sign change in Q{L) then follows from the fact that there 
are more protons than neutrons. In the unfavourable case (bottom), Xtt and Xu are 
interchanged, thus protons favour the prolate shape and neutrons the oblate one in 
high-spin states (the ground state remains prolate). The behaviour of the individual 
contributions follow from the same reasoning. Since there are more protons, they 
dominate the quadrupole moment, and hence the nucleus stays in prolate shape. 

The beneficial effect of ey < in hastening the shape transition can be understood 
in a similar manner. It increases the deformation of protons and reduces that of 
neutrons, leading to a faster change in proton contribution to Q{L) and reducing the 
contribution of neutrons. In the more intuitive picture of Fig. 1 (bottom) , including 
Sy — —0.4 MeV in the Hamiltonian, reduces the energy difference AE between the 
two minima from 0.85 MeV to 0.68 MeV, hence making the cross-over from the 
prolate minimum to the oblate one easier. 

These results for prolate transitional nuclei can be readily extended to oblate 
ones by noting that changing the sign of both Xn and Xu changes the sign of the 
quadrupole moments but their absolute values remain the same. Thus in this case, 
dynamic oblate-prolate shape transitions would be possible when Ny and Xv have 
the opposite signs. The effect of €y on Q{L) remains similar; for A^y > 0, ey < 
pushes the shape transition to an earlier spin, while for A^y < 0, ey > does the 
same. 

We use the above criteria in searching for a potential candidate among the tran- 
sitional nuclei that might exhibit such a dynamic shape transition. The F-spin 
vector parameters in Os-Pt isotopes were recently determined from a study of Ml 
properties, resolving the long standing anomalies observed in the g^-factors of these 
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nuclei Since A^v < in all the isotopes considered in Ref. 0, a negative value 
of Xv in Os isotopes and a positive one in Pt isotopes is a necessary condition for a 
dynamic shape transition. In addition, a positive value of would help it to occur 
at an earlier spin. Inspection of Table 1 in Ref. shows that ^^^Os provides the 
most favourable case for a prolate-oblate shape transition. In Fig. 4, we show the 
level energies, E2 transition m.e. and quadrupole moments in the ground band of 
^^^Os, obtained using the same parameters as in Ref. [^. Although the calculated 
quadrupole moments do not actually change sign, they depict the same rapid change 
seen in Fig. 3, and come very close to doing so. The failure is due to the relatively 
low boson number (A^ = 8), which causes the boson cutoff effect to kick in early. 
The data on E2 transitions (Fig. 4) follow the axial rotor results closely, and do not 
show any sign of boson cutoff. Thus the premature decrease in the E2 m.e. is really 
a problem of the sd-IBM., and for a proper description of high-spin states one needs 
to include (^-bosons in the basis. Unfortunately, an exact sdg-YQM.-2 calculation for 
this nucleus is not feasible at present due to the large basis space. An approximate 
calculation in a truncated basis space is possible but, as a recent study indicates 
||12|| , truncation leads to unreliable results for spins L > 2N. This is exactly where 
the sdg-lBM-2 results are needed, therefore truncated calculations will not serve 
a useful purpose for extending the present results to higher spins. Nevertheless, 
both the proton-neutron deformation difference and the hexadecapole effects drive 
the system in the same direction (from prolate to oblate shape), and inclusion of g 
bosons could only enhance the shape transition. Thus combination of the present 
sd-lBM-2 results with the earher sdg-lBM-1 ones [|], |13], strongly suggests a shape 
transition in the spin range L = 10 to 20 for ^^^Os. The measured quadrupole mo- 
ments deviate markedly from the axial rotor results, but to ascertain whether they 
actually decrease in absolute value would require more accurate measurements. 

In the absence of any other guiding principle, we used the vector parameters that 
were determined from the Ml properties in our search for shape transitions. While 
the F-spin breaking mechanism in the IBM-2 appears to give a mostly consistent 
description of Ml properties in rare-earth nuclei, it is also argued that Ml data 
are really sensitive to single particle degrees of freedom and should not be used to 
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determine collective variables |jT5[ . From this point of view, observation of dynamic 
shape transitions (or lack of it) would help to determine the vector parameters in 
an independent way, and hence it would provide a significant test for the collective 
model interpretation of the Ml observables used in the IBM-2. Structurally, the 
quadrupole parameters used for ^^^Os {xn = —0.68, Xi/ = 0.32) imply competing 
prolate and oblate shapes for the proton and neutron fluids, respectively. This 
picture is similar to the pairing-plus-quadrupole model calculations |T1|, which first 
predicted the static prolate-oblate shape-phase transition in Os-Pt nuclei. 



The errors in the measured quadrupole moments of Os-Pt isotopes ||Tl| are gen- 
erally too large to reach a definite conclusion whether they actually decrease in 
absolute value with increasing spin. More precise measurements of the quadrupole 
moments extending to higher spins are needed to test this prediction. Such measure- 
ments are now possible with the new An high-resolution detector systems "Euroball" 
and "Gamma-Sphere" . It would be very interesting to see whether collective nuclei 
actually exhibit dynamic shape transitions driven by differences between macro- 
scopic shape variables. Apart from that, they would also provide an important 
consistency check on the F-spin breaking mechanism used in explaining the Ml 
properties in the IBM-2. 
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Figure captions 



Fig. 1. Contour plots of the energy surface (8) in the — Pi, plane for a typical 
deformed nucleus (top) and transitional nucleus (bottom) with N^^ = 8, N^, = 4 
bosons. The parameters are (in MeV except for x) = 0, k = 0.15, ^ = 0.17, 
Xv = 2 for both, and Eg = 0.2, Xs = —0.5 in the deformed and Es = 0.4, Xs = —0.2 in 
the transitional case. The prolate (absolute) minima are in the first quadrant and 
the oblate ones in the third. The contour lines are separated by 1 MeV, and the 
absolute minima are -Emm = —4.83 MeV (top), -Emin = —3.07 MeV (bottom). 

Fig. 2. Systematic study of the yrast quadrupole moments in a transitional nucleus 
against the F-spin vector parameters; A^v (a), Xv (b), and e^ (c). The scalar pa- 
rameters are as in Fig. 1 (bottom). The vector parameters are Xv = 2, = (a), 
Ny = 4, Ey = (b), and iVy = 4, Xv = 2 (c). The lines are drawn to guide the eye. 

Fig. 3. Proton-neutron decomposition of the yrast quadrupole moments for the two 
cases with Xv = =t2 in Fig. 2b. The other vector parameters are A^v = 4 and = 0. 
The proton and neutron contributions are indicated by open circles, and the total 
quadrupole moment by filled circles. The lines are drawn to guide the eye. 

Fig. 4. Comparison of the energies, E2 transitions and quadrupole moments of the 
ground band states in ^^^Os calculated in IBM-2 (solid lines) with the experimental 
data (circles) |]11|. The parameters (taken from Ref. [5]) are; N.,^ = 3, N^, = 5, 
Es = 0.4, Ey = 0, K = 0.15, Xs = -0.18, Xv = -1, ^ = 0.17 (in MeV except for 
and x); and e,r = = 0.15 eb. The results for an axial rotor with Qo = 4.6 eb are 
also shown for comparison (dashed line). 
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